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SCATTERING RESONANCES AS VISCOSITY LIMITS 


MACIEJ ZWORSKI 


Abstract. Using the method of complex scaling we show that scattering resonances 
of —A + V, V £ are limits of eigenvalues of —A + V — iex 2 as e —» 0+. 

That justifies a method proposed in computational chemistry and reflects a general 
principle for resonances in other settings. 


1. Introduction and statement of results 

In this note we show that scattering resonances can be defined as viscosity limits, 
that is limits of eigenvalues of Hamiltonians suitably regularized as infinity. The de¬ 
tailed proofs are presented in the simplest case of the Schrodinger operator with a 
compactly supported potential and rely only on standard techniques. 

We consider 

P~~ A + V, Vei“ mp (R"), 

where L“ mp denotes the spaces of bounded functions vanishing outside of some com¬ 
pact set. (Similarly the subscript L* oc denotes function in the space L* on compact 
sets.) The scattering resonances are defined as the poles of the meromorphic continu¬ 
ation of resolvent: 

R v (z) := (—A + V - z)- 1 : L 2 (R n ) -> L 2 (R n ), Imz > 0, 

from the upper half-plane through the continuous spectrum. More precisely, 

R v(z) : L 2 comp (R n ) -A L 2 oc (R n ), (1.1) 

continues meromorphically to the double cover of C when n is odd and to the loga¬ 
rithmic cover of C when n is even. The poles of this continuation coincide with the 
poles of the scattering matrix for the potential V. Their multiplicity (except at the 
threshold 2 = 0) are given by 

771 ( 2 ) := rank j) R v (()d(, (1.2) 

where the integration is over a small circle around 2 - see [DyZ2, Chapter 3]. 

Equivalently, we can consider Green’s function, that is the integral kernel of Ry(z ), 

Rv(z)f(x) = f G(z, x, y)f(y)dy, (1.3) 

Jkp 


1 


2 


MACIEJ ZWORSKI 



FIGURE 1. An illustration of the results of Theorem 1 in the case of 
a specific potential shown on the left. Resonances are computed us¬ 
ing squarepot.m [BiZ]. The eigenvalues of P £ , £ = 1/4 and e = 

1/10 are computing by approximating the unitarily equivalent opera¬ 
tor £^{D 2 + e~^V(e~ 2 x) — ix 2 ) using a basis of eigenfunction of the 
harmonic oscillator D 2 + x 2 . 

and look at the poles of the continuation of z ha G(z, x, y ) for x and y fixed. Another 
way, based on the method of complex scaling, will be reviewed in §2. 

We now consider a regularized operator, 

P £ : = —A + V — iex 2 , £ > 0. (1.4) 

It is easy to see (with details reviewed in §4) that P £ is an unbounded operator on 
L 2 (M n ) with a discrete spectrum. We have 

Theorem 1. Suppose that {zj{£)}ff l are the eigenvalues of P £ . Then, uniformly on 
compact subsets of {z : —7 t/4 < argz < 77 t/4 }, 

Zj(e) —y Zj , £ —y 0T, 

where Zj are the resonances of P. 

Remarks. 1. A more precise statement involving continuity of spectral projections 
is given in §5. The term viscosity is motivated by the viscosity definition of Pollicott- 
Ruelle resonances given in Dyatlov-Zworski [DyZl] - see Example 3 below. 

2. When e < 0 the spectrum of P £ is given by complex conjugates of the spectrum of 
P- e - Hence we have 

Zj (e) —y Zj , £ —y 0—, (1.5) 

uniformly on compact subsets of {z : —7ti/A < arg z < 7T/ 4} . 

3. The term —i£x 2 is an example of a complex absorbing potential and other potentials 
can also be used - see the discussion below. The proof here requires some analyticity 
properties of the complex absorbing potential. 
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4. The restriction to arg z > —7t/ 4 when using —iex 2 is due to the fact that for V = 0 
the spectrum of —A — iex' 2 is given by e^e~’ Kl ^{2\i\ + n), IgN" which is a rescaled 
spectrum of the Davies harmonic oscillator - see §3. One can expand the range using 
ee~ ia x 2 , 0 < a < 7T in which case we recover resonances with aigz > —a/2. 

5. The proof applies with simple modifications to compactly supported black box 
perturbations on M n introduced in [SjZ] - see [DyZ2, Chapter 4] and [Sj3]. In that 
case we need to replace —iex 2 by —ie( 1 — xi x )) x2 where x £ C£°(M n ) is equal to 1 on 
a sufficiently large set - see Example 2 below. 

The computational method based on calculating eigenvalues of P £ was introduced 
in physical chemistry — see Riss-Meyer [RiMe] and Seideman-Miller [SeMi] for the 
original approach and Jagau et al [JZBRK] for some recent developments and refer¬ 
ences. However no rigorous mathematical treatments seem to be available and some 
new interesting open questions can be posed - see Example 4 below. 

Fixed complex absorbing potentials have already been used in mathematical lit¬ 
erature on scattering resonances. Stefanov [St] showed that semiclassical resonances 
close to the real axis can be well approximated using eigenvalues of the Hamilton¬ 
ian modified by a complex absorbing potential. Nonnenmacher-Zworski [NZ1],[NZ2] 
used fixed complex absorbing potentials to study resonance problems employing gluing 
techniques of Datchev-Vasy [DV1],[DV2]. Yet another application was given by Vasy 
in [V] where microlocal complex absorbing potentials were used to obtain Fredholm 
properties and meromorphic continuation of the resolvents (see also [DyZ2, Chapter 
5]). 

We conclude this section with some examples to which Theorem 1 does not apply 
directly but which fit in the same framework. 


Example 1. As explained in [Sj2, (c.31)-(c.33)] the theory of Helffer-Sjostrand [HeSj] 
applies to the case of potentials which are homogeneous of degree m and satisfy the 
condition V(x) =0, i / 0 =>■ Vh(i) ^ 0. That means that resonances of P = 
—A + V can be defined in {z G C, arg z > — 6 * 0 } for some 9 0 > 0. It is interesting to 
ask if the viscosity limit gives a global definition in that case. 

That is easily seen in the case of quadratic potentials. In fact, suppose that 
V ( x ) — + ■ ■ ■ T \ r x r — fi^x r _|_i — ■ • • — ix n _ r x n , A j, > 0. 

As recalled in §3 the eigenvalues of P e , e > 0, are given by 

r n—r 

( A ? - i£ )h 2k j +!) - * Y&j ~ i£ )H 2 k j+r + i), k £ No, 

3 =1 3 =1 
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where the branch of the square root is chosen to be positive on M + . As £ —> 0+ we 
obtain the globally defined set of resonances: 

r n—r 

'y ^ Xj(2kj + 1) — i ~y Hj(2kj +r + 1), k G Ng. 

3 =1 J=1 

Example 2. This example fits in the framework of black box scattering with one 
dimensional infinity. Consider the modular surface M = SX 2 (Z)\HI 2 and A m A 0 
the Laplacian on M. We then put P = —A m ~ \ where \ guarantees that the the 
continuous spectrum of P is given [0, oo). This is a black box Hamiltonian on "Ho © 
L 2 ([0, oo)) in the sense of [SjZ] - see [DyZ2, §4.1], Traditionally, the resonances of the 
quotient M are defined as poles of the meromorphic continuation of (-Am-s(I-s))^ 1 
from Re s > \ to C and are given by the embedded eigenvalues with Re s = 4 and by 
the non-trivial zeros of C(2s) where ( is the Riemann zeta function. The resonances of 
P are then given by s(l — s). 

If we choose the fundamental domain of SX 2 (Z) to be {x + iy : \x\ < l,x 2 + y 2 > 1} 
then the Laplacian in the cusp y > 1 is y~ 2 (d 2 + d 2 ). The Hamiltonian on L 2 ([0, oo) r ) 
is given by — <9 2 , r = log y - see [DyZ2, §4.1, Example 3]. In the language of Theorem 1 
(see Remark 5) and in (x, y) coordinates 

P £ = -A M ie( 1 - x{y))(\°gy) 2 , 

where x G C^°([0,oo)), y(y) = 1 for y < | and y(y) = 1 for y > 2. The operator P £ 
has discrete spectrum for £ > 0 and the eigenvalues converge to the resonances of P 
uniformly on compact subsets of Irn z > —7 t/4 . Equivalently if we define E £ 

s(e) G E e -«=> s(e)(l - s(e)) G a(P e ) 

The limit points of E e , £ —>■ 0+, in Res < |s| > C are given by the nontrivial zeros 

of C(2s). 

Example 3. Suppose that A" is a compact manifold and V is a vector field on M 
generating an Anosov flow, ^ = exp tV. That means that the tangent space to A" 
has a continuous decomposition T x X = E 0 (x) © E s (x ) © E u (x) which is invariant, 
dtp t (x)E,(x) = E,((f t {x)), E 0 (x) = iy(x), and for some C and 9 > 0 fixed 

<Ce~ m \v\ x , v G E u (x), t < 0, 

< Ce~ ew \v\ x , v G E s (x ), t > 0. 

where | • \ y is given by a smooth Riemannian metric on X. A class of examples is 
given by X = T l Ad where M is a negatively curved Riemannian manifold and yA is 
the geodesic flow in its unit tangent bundle X. 
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If A g < 0 is the Laplacian for some metric on X then - see [DyZl] - the limit set of 
the spectrum of 

P £ — X/i -f- isXg, £ —0+ 

is a discrete set given by the Pollicott-Ruelle resonances - see [DyZl] for the defini¬ 
tion and references. Adding the Laplacian corresponds to taking a viscosity regular¬ 
ization and that explains our terminology. Another interpretation is given in terms 
of Brownian motion: the pullback by the flow flow x(t) := <p f (o;(0)), is given by 
e ltp °f(x) = f(x(t)), x{t ) = #(0) = x. For e > 0 the evolution equation is 

replaced by the Langevin equation: 

e~ UPe f(x) = E [/(z(t))], x(t) = -14(f) + y/2eB(t), x(0) = x, 

where B(t) is the Brownian motion corresponding to the metric g on A". Hence con¬ 
sidering P £ corresponds to a stochastic perturbation of the deterministic flow. In the 
case of scattering resonances the same interpretation can be proposed on the Fourier 
transform side. 

The assumption that the flow satisfies (1.6) is crucial as otherwise the limit set is 
typically not discrete. The simplest example is given by A" = S 1 x S 1 , S 1 = M/27tZ, 
and V = d Xl + ad X2 , a ^ Q, A g = d x + <9 2 2 . In that case the limit set of the spectrum 
of P £ is the lower half plane. Other limit sets are possible, for instance in the case of 
the geodesic flow on § 2 , X = T 1 § 2 ~ SO( 3). The spectrum of Po is given by Z (with 
infinite multiplicities) and if we take A g to be the Casimir operator then the limit set 
of the spectrum of P £ as e —>■ 0+ is Z — i[0, oo). For yet another example see [DyZl, 
§!]• 


Example 4. We expect that viscosity definition of resonances remains valid, in a small 
angle near the real axis, for all dilation analytic potentials - see [HeSj] and references 
given there and §2 below for a review of complex scaling. It would be interesting to 
find a Schrodinger operator P for which the limit set of the spectrum of P £1 £ —>• 0 is 
not discrete. Candidates are given by potentials which are not dilation analytic, for 
instance, 



sin x 

1 

X 




Notation. We use the following notation: / = On(g)H means that ||/||# < Cgg where 
the norm (or any seminorm) is in the space H, and the constant Ca depends on t. 
When either i or H are absent then the constant is universal or the estimate is scalar, 
respectively. When G = Ge(g) : H\ — > H 2 then the operator G : Hi —> H 2 has its norm 
bounded by C^g. Also when no confusion is likely to result, we denote the operator 
/ t-A gf where g is a function by g. 



6 


MACIEJ ZWORSKI 


Acknowledgments. The author would like to thank Mike Christ, Semyon Dyatlov, 
Jeff Galkowski, John Strain and Joe Viola for helpful discussions. This project was 
supported in part by the National Science Foundation grant DMS-1201417. 


2. Review of complex scaling 


The complex scaling method changes the original Hamiltonian P = P 0 to a non-self- 
adjoint Hamiltonian P o ,0 such that P o ,0 — 2 : H 2 At 2 is a Fredholm operator when 
arg z > — 26 . It was introduced by Aguilar-Combes [AgCo], Balslev-Combes [BaCo] 
and Simon [Si]. For a review of practical applications of this method in computational 
chemistry see Reinhardt [Rci], As the method of perfectly matched layers (PML) it has 
reappeared in numerical analysis - see Berenger [Be] . The presentation here follows the 
geometric approach of Sjostrand-Zworski [SjZ], Eventually the proof that the viscosity 
eigenvalues converge to scattering resonances is a straightforward application of the 
methods of [SjZ] (see also [Sj3, §7.2] for a more detailed presentation and [DyZ2, §4.5] 
for an approach to complex scaling based on the continuation of the Green function 
G(z,x,y ) in (1.3) in variables x and y). 

Suppose that O C C" is an open subset and that 

P(A D z ) = a a(z)Df, D z . : = \d z ., Df = Dfl ■ ■ ■ (2.1) 

|a|<m 

is a differential operator with holomorphic coefficients. For instance we can have 
P(z,D z ) = Z% 1 Dl-tez]. 

Suppose that fl C €" is an open subset and that T C 11 is a maximal totally real 
submanifold. That means that T is a smooth real submanifold of dimension n such 
that 

Vi 6 T, T x rniT x T = {0}. (2.2) 

Here we identify T X Y with a real subspace of C n . The condition (2.2) means that there 
exists a complex linear change of variables A : C n —>■ C" such that A(T X T) = M n C C n . 
Locally, V can be represented using real coordinates: 


R n DU3x^ f{x) = ifiix),--- , f n (x)) ercHc C n . (2.3) 


Composing the matrix d x f(x) := (d Xj fk{x)) i<k,j<n with A we obtain an invertible 
matrix M n —y M n . That means that 


( dfkjx) 

V 9xj 


t^O. 


det 


l<k,j<n 


(2.4) 
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Conversely, if (2.4) holds, then df(x) is an injective complex linear matrix and for any 
sets U, V C C n , df(x)(U) fl df(x)(V) = df(x)(U D V). Hence, 

T X T n iT x r = df(x)(R n ) n idf(x)(R n ) = df(x)(R n ) n df(x)(iR n ) 

= df(x)(R n niR n ) = {0}, 

and (2.4) implies (2.2). The volume form on T is obtained by pushing forward the 
standard volume form on M n by /. That of course depends on the choice of / (in what 
follows the uniformity will be guaranteed by (2.8) below). 


Example. 

Then 


As a simple illustration consider n = 2 and f(x i, x 2 ) = {x\ + ix 2 , 0) G C 2 . 


d x f(x) 


T x f(R 2 ) = C®{0} C C 2 . 


The tangent space is not totally real and condition (2.4) is violated. To introduce the 
next topic we also note we cannot restrict operators, P, with holomorphic coefficients 
to f(R 2 ) in a way that for holomorphic functions, u, (Pu)|/(r 2 ) = (P/(r«))(m|/( irp)- As 
an example consider P = d Z2 and u = z 2 . 


The point of introducing totally real submanifolds T is the fact that an operator, P, 
with holomorphic coefficients can be restricted to an operator with complex smooth 
coefficients on T, Pp, in such a way that for u holomorphic near T, Pu\p = -Pr(w|r)- 

The differential operator P(z,D z ) given in (2.1) defines a unique Pp a differential 
operator on T as follows. Using (2.3) we can identify a small neighbourhood of any 
z 0 G T with U C R n . Then u G C°°(T fl f(U)) can be identified with mo/g C°°(U). 
We then have 

(Ppu) o f( x ) = (a a o f)(x)(( t d x f(x)~ l D x ) a (u o f)(x). (2.5) 

|a|<m 

It is easy to see that this definition is independent of the choice of / and that the 
condition (2.4) is crucial. 

The key fact is the standard result about continuation of solutions to Ppu. The 
proof based on [Le], [Ma] and [Sjl] can be found in [SjZ, Lemma 3.1] and (in more 
detail) [Sj3, Lemma 7.2], With the notation above we have the following: 

Lemma 1. Suppose that W C M n is open and that F : [0,1] x W 3 (s,x) i-G F(s,x ) G 
C n , is a smooth proper map satisfying for all s G [0,1] 

det d x F(s, x) ^ 0, and x i—>• F(s, x) is injective. 

In addition assume that there exists a compact set K C W such that 

x eW\K ==>• F(0, x) = F(s, x), 0 < s < 1, 
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and that F([0,1] x W ) C with P(z,D z ) a differential operator with holomorphic 
coefficients in fh 

Now assume that for T s := F({s} x W), Pr s is an elliptic differential operator in 
the sense that 

| X > C|cr, (z,()err,. 

\a\=m 

If uq £ C°°(r 0 ) and Pr 0 uo extends to a holomorphic function on Q, then for every 
s £ [0,1] there exists a holomorphic function, U s defined near such that, for some e, 

Uo\r 0 = u 0 , |s — s' | < £ ==>- U s = U s ’ on the intersection of their domains. 

In other words, the function uq defined on r 0 extends to a possibly multivalued 
function U in a neighbourhood of /([0,1] x W). 

The lemma will be applied to a family of deformations of in C n . Our goal is to 
restrict the operator P £ = —A — iex 2 + V, £ > 0, to the corresponding totally real 
submanifolds. For that the deformation has to avoid the support of V and we choose 


r 0 such that suppO C fi(0, r 0 ). We then construct 

[0,7r) x [0, oo) 3 (9,t) i —> g B (t) e C (2.6) 

which is C°°, is injective on [0, oo) for every fixed 6 and satisfies 

go(t) = t for 0 < t < To, (2.7) 

0 < arg g 0 (t)<9, d t g 0 0, (2.8) 

ar g ge(t) < arg d t g 0 (t) < arg g e (t) + e 0 , (2.9) 

ge(t ) = e ie t for t > T 0 where T 0 depends only on e 0 and r 0 . (2.10) 


We now define the totally real submanifolds, Tq, as images of M n under the maps 
fo-.R n ^C n , f e (x):=g e (\x\)x/\x\, T e := f e (R n ). 

For £ > 0 and 0 < 6 < n we put 

-A e := (A 2 )|r e , x e := z\r e , 

Qs,e ■= -A g - iex 2 e , P efi := Q £ ,e + Id- 

Parametrizing r,? by (t,ui) £ [0, oo) x § n_1 , (t,u) ha- gg(t)u>, we have 

-A 0 = (g' e {ty l Df) 2 - i(n - 1 )gg(t)~ 1 g'g(t)~ 1 D t + g e (t)~ 2 Dl , 
where D t = dt/i and D‘f = — A§n-i. The symbol is given by 

o’(-Afl) = g' e (t)~ l T 2 + g e {ty 2 w 2 , ( t,u;r,w ) £ T*([0,oo) x § n_1 ). 

The basic result based on ellipticity at infinity is 

-26+6 < mgz < 2w-20-S, \z\ > 6 =+ (-A e -z)-' = 0,(1) : L 2 (r«) -> V 2 (r s ). 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 
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This follows from [SjZ, Lemmas 3.2-3.5] applied with P = —A. As will be reviewed 
in §4 this shows that P 0 ,6> — z : H 2 —> L 2 is a Fredholm operator in this range of values 
of z and that the eigenvalues are independent of 6. 

The crucial property is 

Lemma 2. Let R 0 (z) = (—A — z )^ 1 : L 2 —» H 2 , Irn z > 0, be the free resolvent and let 
Ro(z) also denote its analytic continuation across [0,oo) as an operator L 2 omp —» H 2 oc . 

Suppose that x G C^°(P(0,r 0 )) so that x is defined on Tg. Then for —26 < arg z < 
2tt -26, 6 < tt, 

xRo(z)x = X(~^e ~ zY l x- (2.14) 

Proof. We recall the main features of the proof which is implicit in [SjZ, §3]. It is 
sufficient to establish the identity (2.14) for 0 < argz <2ix — 26 as it then follows by 
analytic continuation. It is also enough to show that in this range of z and 0 < 6\ < 

6 2 <6, |0i-0 2 | <1, 

y(-A 01 - z)~ 1 x = X(-Ae 2 - zY l x■ (2.15) 

For that we show that for / e L 2 (P(0, ro)) C L 2 (Tg j ) there exists U holomorphic in a 
neighbourhood ^e 1 ,e 2 of 

U (r«\fl(o,m))<=C” 

0i<0<0 2 

such that 

U\r g .(x) = [("A 0 . - z)~ 1 xf](x) for x G Tg. \ B(0,r Q ). (2.16) 

The unique continuation property for second order elliptic operators then shows that 

x{A 0l - z)~ 1 xf = x(^e 2 ~ z)~ 1 xf , 

proving (2.14). 

To show the existence of U such that (2.16) holds we use Lemma 1 applied to a 
modified family of deformations. The key is to show that a holomorphic extension, 
I/, of the solution to (—A^ — z)u\ = %/, U\ G L 2 (Tg 1 ), restricts to u 2 G L 2 (Tg 2 ) 
(the equation (—A g 2 — z)u 2 = xf is automatically satisfied). That means that u 2 = 
(-A g 2 - z)- 1 ^f) proving (2.16). 

The modified family of contours is obtained as follows. Fix T iX> l and choose 
X G C£°((2, 5); [0,1]) equal to 1 near [3,4], Then define 

ge 1 } 02 ,T(t) := gefit) + x(t/T)(gg 2 (t) - gefit)), 
re»i,e» 2 ,T : = { 9 Oi, 02 ,T(t)uj : t G [0, oo), u) G § n : } C C n . 

We can apply Lemma 1 to the family of totally real submanifolds interpolating be¬ 
tween T 0l and T 0 u e 2 ,T- [0,1] 3 s 1 —> ^e 1 ,ei+s{e 1 ,e 2 ),T- That implies that there exists a 
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FIGURE 2. The deformed totally real submanifold Fo, q 2: t interpolating 
between T^ and Ve 2 . 

holomorphic function U T defined in a neighbourhood of the union of these submani¬ 
folds and such that u\ = U T |r 9 . Changing T we obtain a family of functions agreeing 
on the intersections of their domains and that gives U dehned in the neighbourhood 
YLq 1 q 2 . To see that U\t ( , 2 G L 2 {Y e2 ) it suffices to show that 

\\U T \re 1 ,g 2:T \\L 2 (rg li e 2 , T ) ^ C'o||«l|U 2 (r ei n{T<|z|<6T},)) ( 2 - 17 ) 

where C 0 is independent of T. (We apply (2.17) with T = 2 J and sum over j.) 

To see (2.17) 

Qi(T) = {z G C n : 2T < |z| < 5T} fl Tq 1 0 2 T D Y dl ^ 2)T \ r 01 , 
n 2 (T) = {z EC n :T <\z\< 6 T} n T 6 u e 2 , T , f) 2 (T) \ Oi (T) C e i6l M n . 

We claim that for T large and u G C°°(re lt e 2 ,T), 

||w||z, 2 (f2l ( T )) < C'lK—Ar flii g 2iT - z)u\\ L 2(n 2 (T)) + C'll' u IU 2 (n 2 (T)\n 1 (T))- (2.18) 

For 102 — | 1, this estimate is a perturbation of a standard semiclassical ellip¬ 

tic estimate: treating h := 1/T as a semiclassical parameter, uniform cllipticity of 
—e~ 2 ,e h 2 A — z shows that for v G C°°(M n ), 

IMU 2 ({2<|a:|<5}) A c || (—e 2l9 lrA — ^)v||z, 2 ({ 1 <|a:|<6}) + C||u|| L 2 ({ 1 <| s |< 2 }. U { 5 <|. E [< 6 }. 
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(This can be seen applying the inverse from [Z, Theorem 4.29] to x v where x £ 
Q°((l, 6)) is equal to 1 on [2, 5].) The properties of Llj(T) then imply (2.17) completing 
the argument. □ 


3. The Davies harmonic oscillator 

The operator H ea := —A + e _ * 7 £x 2 , £ > 0, 0 < 7 < tt, was used by Davies [Da] to 
illustrate properties of non-normal differential operators. We recall the following basic 
result: 

Lemma 3. The operator H en is an unbounded operator on L 2 with the discrete spec¬ 
trum given by 

a{H £}1 ) = e- i 'i /2 V^(n + 2\N™\), |a| = a x + • • • + a n . (3.1) 

If Q <1 {z : —7 < argz < 0} \ e _ * 7 / 2 [0, oo ) 7 then for some constant C\ = C'i(fi), 

^~ h/Cl <\,ze a (3.2) 

C'l 

In addition for any S > 0 there exists a constant C 2 such that, uniformly in e > 0, 
\\{H e ^ - zY l \\ L 2^ L 2 <c 2 /\z\, 8 < argz < 27r - 7 - <5, \z\> 5. (3.3) 

Proof. By rescaling y = y/ex this operator in unitarily equivalent to —cA y -\-e~ ll y 2 , that 
is a semiclassical, h = yje, quadratic operator. For the analysis of the spectrum and 
upper bounds on the resolvent for general quadratic operators see Hitrik-Sjostrand- 
Viola [HSV] and references given there - in particular we obtain (3.1) and the upper 
bound in (3.2). The lower bound in (3.2) follows from general arguments for operators 
with analytic coefficients - see [DeSZ, §3] and the bound (3.3) from (semiclassical) 
cllipticity of —h 2 A y + e~ ll y 2 — z for 8 < arg z < 2n — 7 — |z| > 8. □ 

We now consider the special case of H s ^/2 = Qe, 0 and of its deformation Q £) e - see 
(2.12). The facts we need are given in the next two lemmas. The first is the analogue 
of Lemma 2: 

Lemma 4. In the notation of Lemma 2, 0 < 9 < 7 t/ 8 , e > 0, and —29 < arg 0 < 
37t/2 + 26 we have 

x(Qe ,0 - z)~ l X = x(Qs,e ~ z)~ x x- (3.4) 

In particular, for 0 < 6 < 7 r /8 , the spectrum is independent of 6 and given by 
yfee-^^fnP 2|Ng|). 

Proof. We follow the argument in the proof of Lemma 2 and use the notation intro¬ 
duced there. Hence it is enough to prove that 0 < 9i < 9 2 < tt/8 and \9i — 9 2 \ small it 
is enough to show that 


x{Qe,e 1 - z) l x = x(Qs,e 2 ~ z) l x ■ 
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Figure 3. A visualization of the spectrum of Q £) q = —A — lex 1 which is 
equal to the spectrum of the deformed operator Q £} e- The lightly shaded 
region is the numerical range of Q £y o and the darker shaded region, the 
numerical range of — e~ 2l6 A — ie 2 l 6 ex 2 . The estimates for the resolvents 
°f Qe,e improve outside of that region. 

We only need to establish this for z e e n- 2 0i+ 7r / 2 )(i ) qq) as then the result follows by 
analytic continuation. The only difference is an estimate which replaces (2.18): for 
r > 1, 

IMU 2 (fii(T)) < c \\(Qr g ^e 2 ,T ~ ^ e ” 201r )^IU2(f2 2 (T)) + C'||«||L2(a 2 (T)\n 1 (T)), ^ 

Qdi, 62 ,T ■ g lt e 2 ,T g 1 : g 2tT ) 

uniformly for T 1. To see this we first note that for e > 0, Qo 1 ,b 2 ,t — z, z G C, is a 
Fredholm operator (since it is elliptic and near infinity it is equal to e~ 2 l 0 H £j 7 r / 2 - 4 . 9 ) ■ 

To obtain an estimate we notice that for t > T 

(t) = X(t/T)e^ + (1 - x(t/T))e + (t/T) x '(t/T)(e M * - e ife ), 
so that from ( 2 . 8 ) and ( 2 . 10 ), 

9\ - C\0 2 - 0i| < arg g' 6 ueaiT (t) < 9 2 . 

Also, d 1 < arg ge u 0 2 , T (t) < 9 2 . Hence, 

Re((e 2 ldl Q dl ,e 2 ,t ~ ir)u,u) > \\Du\\ 2 /C 
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where we used the fact that for for 0 < 6 < 7 t/ 8 , Re(— ie Ae ) > 0. The imaginary part 
is then estimated as follows, 

-Im ((e 2l9l Qo lt o 2 , T ~ ir)u,u) > T\\u\\^r ei ,e a ,T) - 0(|0 2 - 9i\)\\Du\\ 2 . 

We conclude that when |0 2 — 0i| is small enough 

II«IW ~ ie- 2ie <T)u\l > (|| U || + ||Du||)/C, 

This and the Fredholm property imply that 

(Qe i,0 2 ,t — 2,9l T ) 1 = 0(1) : L 2 (T 6li e 2>T ) —> H 1 (Tg lj92tT ). 

that is the operator is invertible with bounds independent of T. From this (3.5) 
follows by a standard localization argument: we choose Xt € C°°(Q, 2 (T), [ 0 , 1 ]), such 
that xt — 1 on fli(T) with derivative bounds independent of T. We then apply the 
inverse above to (Qe 1 ,e 2 ,T — ie~ 2l6l T)xTU with the commutator terms estimated by 
IM|L2(Q 2 (T)\ni(T))- □ 

The next lemma shows how complex scaling dramatically improves the exponential 
bound (3.2): 

Lemma 5. Suppose that 0 < 6 < 7 t /8 and that Q d {z : —26 < arg z < 37 t /2 + 26}. 
Then there exists C = C(f2) (in particular independent of e > 0) such that 

|| (Qe,d ~ z) 1 1|X/2—>Z,2 < C, ZGfl. 

Proof. Let Xj G C£°([0, oo)) be equal to 1 on [0, r 0 ] and satisfy Xj — 1 on supp Xj+i , 
j = 0,1. Parametrizing Tg by Fq : [0, oo)t x S n_1 —* Tg, Fg(t,u >) = go(t)u (with go 
given in ( 2 . 6 )) we define functions Xj £ as 

Xj° F 0 (t,u) :=Xj(th), 0 < h < 1. 

In view of (2.10) and (2.13) we see that for h small enough 

Q m (1 - xf) = (-e^A, - iee™x 2 )( 1 - *f) 

= e~ 2l6 H £n (l~Xi), 7 := vr/2 — 40, x = tu. 

In view of (3.3) we have 

(1 - x h 2 )e 2ie (H £n - e 2ie z)-\l - X h 2 ) = O s { 1) : L 2 (V e ) —>• H 2 (T e ), (3.6) 

for 

— 5 < 26 + arg z < 2 tt — 7 — 5 = 37 t/2 + 40 — <5, |z| > 5, 

and in particular for z 6 fh We stress that the bounds are independent of e. 

Noting that 


(— A» - z )- 1 = 0(1) : i 2 (r») -> // 2 (U), 2 e SI, 


( 3 . 7 ) 
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(for 0 < 9 < 7t/8, —29 < arg z <2n — 29) we now put 

T e h e (z) := X h 0 (-A e - z)- 1 ** + (1 - x\)e 2i \H £ „ - e 2ie z)~\l - **), 
so that (Q s ,0 — z )T £ g(z) = I + K^q(z), where 

K s,ei z ) '■= - iexlxo(-Ao ~ ^) _1 Xi - [A*,Xo](-A 0 - z )-\i 

+ [A 0 , xfc^l - x*)(tf £ , 7 - - Xa)- 

Since [A fl , X J] = O(h) : H\Y e ) -A L 2 (r 0 ) and x 2 X ? = ^(^ 2 ) : F 2 (r,) —> L 2 (r e ), we 
conclude from (3.6) and (3.7) that for z G fi, 

K%(z) = 0(h~ 2 e) + 0(h) : L 2 (r 0 ) — >■ L 2 (r 0 ). 

Hence by choosing h first, we see that for e < e 0 (h), I + K^ e (z) has a uniformly 
bounded inverse and 0 < e < £ 0 

{Qe, h - z)- 1 = T^ g (z)(I + K^(z))- 1 = 0(1) : L 2 (r 0 ) —>■ L 2 (T e ), z G Q. 

In view of Lemma 4 we know that for z G O, (Q £) h — z)~ l exists for e > e® and that 
gives the bound for all values e. □ 


4. Meromorphic continuation 

In this section we will review the meromorphy of the resolvent Ry(z), see (1.1), in 
a way connecting it to the resolvent of P £ given in (1.4), e > 0. For that we define 

R e (z) = (—A — iex 2 — z) -1 , Rv, e (z) = (—A — iex 2 + V — ^) _1 , £ > 0. (4.1) 

For e > 0, these operators are meromorphic for z G C as operators on L 2 . For £ = 0, 
Ro(z) is holomorphic in the sense of (1.1) on the double cover of C \ {0} when n is 
odd and on the logarithmic cover when n is even - see for instance [DyZ2, §3.1]. We 
are only concerned with continuation to argz > —7t/4. 

Let p G C£°(R n ; [0,1]) be equal to 1 on a neighbourhood of supp V. We have 

Lemma 6. For £ > 0 

zg (/ + VR £ (z)p)^ 1 , —7t/4 < arg 2 < 77t/4, 

is a meromorphic family of operators on L 2 (M n ) for with poles of finite rank. Then 

m £ (z) := tr j)(I + VR £ (w)p)~ 1 d w (VR £ (w)p)dw , (4.2) 

where the integral is over a positively oriented circle enclosing z and containing no 
poles other than possibly z, satisfies 
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where m(z) is the multiplicity of the resonance z given by (1.2). 


Proof. We recall the standard argument (see [DyZ2, §2.2, 3.2] and references given 
there). For any 6 > 0 and uniformly in e > 0, 

R e {z) = 0 6 {l/\z\) : L 2 (W) -> L 2 (R n ), S < arg 2 < 3tt/2 - <5, \z\ > 5. (4.4) 


This follows from self-adjointness for £ = 0 and from (3.3) for £ > 0. 
For z in (4.4) and Q e := —A — iex 2 , 

{P e -z) = (Q e -z){I + R e {z)V) 

= (I + VR e {z)p){I + VR £ {z){ 1 - p))(Q e - z 


Noting that 


(/ + VR e {z){ 1 - p))- 1 = I - VR e (z)( 1 - p) 
we obtain from (4.4) and (4.5) that 

RvM = Re{z){I + VR £ {z)p)-\I - VR s (z)( 1 - p)), 
5 < argz < 37t/2 — 5 , \z\ 1, 


(4.5) 


(4.6) 


where for large \z\, I + VR £ (z)p is invertible by a Neumann series argument. Since 
z i—^ VR £ (z)p is a holomorphic family of compact operators for —7 t/ 4 < arg z < 37t/ 4 
(see Lemma 3 for the case £ > 0), z i-> (/ + VR £ (z)p)~ 1 is a meromorphic family 
operators in the same range of z. (For £ > 0 the meromorphy is in fact valid for 
z G C - see [DyZ2, §C.4].) The formula (4.6) remains valid for that range of z with 
boundedness on L 2 for £ > 0. For £ = 0 we note that 


(1 - WJ„M(1 - p)) , (/ + VR„(z)Y' : L 2 comp -4 i c 2 omp , 


floM : r 


2 

comp 


-A L 


2 

loci 


and we obtain the meromorphic continuation of Rvp(z) : L 2 omp —>• L 2 oc . Arguing as in 
the proof of [DyZ2, Theorem 3.23] we obtain the multiplicity formula (4.3). (This can 
also seen using complex scaling as reviewed in the proof of Theorem 2 below.) □ 


5. Proof of convergence 

The proof of convergence is based on Lemma 6 and on the following lemma in which 
we use the complex variable techniques of §§2,3. 

Lemma 7. For x £ C£°(M n ) consider 

T*(z) x(—A ^ iex 2 — z) -1 £ 2 (— A — z) _1 x, 0 < arg z < 3ir/2. (5.1) 

Then Tf continues to a holomorphic family of operators 

T£(z) : L 2 — » L 2 , — 7t/4 < arg z < 77t/4. 
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If D, <1 {z : —7t/ 4 < arg z < 37r/2} then there exists C = Cq >x (independent of e) such 
that 

\\T?(z )\\ L ^ L 2 <C, zen, £>0. (5.2) 

Proof. In the notation of (4.1) we see that for 5 < arg z < 37t/2 — 5, \z\ > 5, 

x(Re(z) ~ Ro(z))x = i£XRs(z)x 2 R 0 (z)x , 

where we note that, for in our range of Rq{z)x '■ L 2 —» e~ cs ' x ^L 2 (by looking, for 
instance at the explicit formulas for the resolvent, see [DyZ2, §3.1], or by conjugation 
with exponential weights) and consequently x 2 Rq(z)x : L 2 —>■ L 2 . Hence 

T f(z) = - l -(xRe(z)x ~ xRo(z)x). (5.3) 

The right hand side is holomorphic for —7 t/ 4 < arg z < 57t/ 4 which provides holomor- 
phic continuation of Tf(z), £ > 0. 

We now use Lemmas 2 and 4. For that we choose r 0 in the definition of large 
enough so that suppy C B( 0, r 0 ) and take 9 = 7t/ 8. Then we have 

T H Z ) = -- £ (x(Qs,e ~ -) _1 X - x(Qo,e ~ -) _ 1 x) ^ 

= x(Qe,e ~ ^) _ 1 ^(Q O ,0 - z)~ l x, 

where, in the notation of (2.12), xe '■= x|r e - We now note that for zed, 

(Qo,e ~ z)~ l X ■■ L\ T 0 ) —>■ e -nl-l L 2 (r e ). (5.5) 

This can be seen by conjugation by exponential weights or by constructing a parametrix 
for Qo ,0 as in the proof of Lemma 4 and using the explicit properties of (— e~ 2l9 A — 
z)^ 1 = e 2%e Ro(e 2ie z) . From this and Lemma 4 we obtain 

II (Qe,e - z)~ 1 x 2 e {QQ )d ~ z)~ 1 x\\l^l^ < Cq, z e fi. 

Inserting this into (5.4) concludes the proof. □ 

We can now state a stronger version of Theorem 1 formulated using the projections 
appearing in (4.2): 

Theorem 2. Suppose that —7t/ 4 < arg 2 < 57t/4 and that m(z) = m > 0, where m(z) 
is the multiplicity of the resonance of P := —A + V at z - see (1.2). 

Then there exists £q and 5 such that for 0 < £ < £ 0 , P £ = —A + V — i£x 2 has m 
eigenvalues in D(z,5): 

trn s = m, n £ [ (( - P^dC, U 2 £ = n e , (5.6) 

zm JdD(z,S) 

and for any x G Q°(M n ), 

X^sX e C°°([0,e 0 ], T(L 2 (M n ),L 2 (M n ))). 


(5.7) 
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Remarks. 1. Notation / G C l {[a,b]) means that /, and f continuous in [a, &]; here 
) are the left and right derivatives of / at those points. By induction we 
then define C k ([a,b]) and C°°([a,b\). In view of (1.5) we cannnot expect analytic 
dependence on e. 

2. For x = 1 on supp V, m(z) = rank x^oX and (5.7) shows the convergence of resonant 
states in the case of simple resonances. As the proof shows a stronger statement is 
obtained by using the complex scaled operators: for 9 = tt/8, 

n £> e := [ (C - Pe^dC, P £ ,e = -A|r e - ie(x |rj 2 + V, 

JdD(z,8) (5.8) 

U £t e eC([0,e o )-C\L 2 (r e ),L 2 (r e ))), n £ #x e C°°([0, e 0 ); /^(L 2 ^*), L 2 (r e ))), 

where Tg is the deformation dehned in (2.11). 

Proof. We first note that (4.6) and Lemma 4 imply that for —7 t/4 < —29 < arg z < 
27r — 2 9, e > 0, 

(Pe,o - z)- 1 = 0 Qs,e - z)-\l + VR e (z)p)-\l - V(Q £! g - z)~\ 1 - p)). (5.9) 

Since z t-A (Q e ,o ~ z)^ 1 is a holomorphic family in our range of z’s, the Gohberg-Sigal 
theory - see [DyZ2, §C.4] - shows that the poles of (P e ,o — z) _1 with arg z > —29 are 
independent ofO<0<7r/8 and 

tr j(P,.9 ~ O' 1 * = tr O-VC. e > 0. 

If in the dehnition of we take r 0 large enough so that supp x C B( 0, r 0 ) then Lemmas 
2 and 4 show that yLL 0 X = xJ^eX- 

Hence it is enough to prove (5.8). For that we note that in the notation of Lemma 7, 

(/ + VR £ {z)p)~ l - (/ + VRo(z)p)- 1 = ie{I + VR £ {z)p)- l Tf{z)(I + VR 0 {z)p)~ l 

= O z (e ||(/ + VR £ (z)p)- 1 ll L2 ^ L2 ) : L 2 -> L 2 . 

We can now apply the Gohberg-Sigal-Rouche theorem [DyZ2, Theorem C.9] to see 
that the poles of (/ + VRo(z)p) and (/ + VR £ (z)p ) _1 coincide with multiplicities. 
This and (5.9) prove the first statement in (5.8). The second statement follows from 
differentiation and estimates similar to (5.5). □ 
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